Static Hopfion solutions of the extended 
Skyrme-Faddeev model 



David Foster* 

Instituto de Fisica de Sao Carlos; 
IFSC/USP, Universidade de Sao Paulo - USP, 
Caixa Postal 369, CEP 13560-970, 
Sao Carlos-SP, Brazil 



Abstract 

We construct non-axially symmetric static soliton solutions, with non-zero topo- 
logical charges, of an extension of the Skyrme-Faddeev model. The model has an 
extra quartic-derivative term and we choose its coupling to the Skyrme-term to 
be negative. We solve the full equations of motion to find numerical solutions 
with topological charge up to seven and find that the model favours large ring-like 
solutions. 



1 Introduction 

The Skyrme-Faddeev model [I] is a non-linear model in three-dimensional space where 
the field takes its value on the unit two-sphere and the minimum-energy soliton solutions 
are called Hopfions. Substantial numerical work, j2j [3] etc, has shown that the model has 
vortex-string like soliton solutions and can resemble either rings, links or knots. Using 
the Cho-Faddeev-Niemi-Shabanov decomposition, Faddeev and Niemi conjectured that 
the Skyrme-Faddeev model describes the low-energy limit of the pure SU (2) Yang-Mills 
theory [4]. This conjecture is somewhat controversial and there is evidence either way, 
also Faddeev has recently proposed a modification [5] . 

When Gies [6] calculated the Willsonian effective-action, to one-loop, for the SU (2) Yang- 
Mills theory he discovered the Skyrme-Faddeev model with an extra quartic-derivative 



term. This model is known as the extended-Skyrme-Faddeev model (2.1), and is believed 
to have a stronger basis in high-energy physics than the Skyrme-Faddeev model. 
The extended-Skyrme-Faddeev model was first investigated in [7J, where using an axially 
symmetric ansatz they found a number of non-trivial Hopfions. This was for a differ- 
ent coupling region to the one investigated here. Later there was another study of the 
extended-Skyrme-Faddeev model [8], here they noticed a different region of stability and 
using an axially symmetric ansatz they found a number of non-trivial solutions. In [8] the 
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authors found a noticeably different topological charge four solution than the one found 
in the Skyrme-Faddeev model. Hence this model requires further study. 
This paper is separated into three sections. In section [2] we shall introduce the extended- 
Skyrme-Faddeev model and discuss the two coupling constant regions, in section [3] we 
briefly introduce a technique [2j which is used to create topologically non-trivial knot- 
ted/linked initial configurations, in section [I] we discuss our charge specific numerical 
solutions for two coupling regimes. Then there are some concluding remarks. 



2 The extended- Skyrme-Faddeev model 

Here we are only interested in static-minimum-energy solutions, so the model is best 
defined by the static-energy functional 

E = ^ (di<f> ■ d i( f> + ^FijFi - ^(d i( j> ■ <9^) 2 ) d 3 x, (2.1) 

where F^ = Cabc^adi^djcj) . The field is a three-component unit vector = (0i, 02,03) 
and is the map : M 3 — » S 2 . Finite energy requires that the field tends to a constant at 
spatial infinity, and we choose this to be 0°° = (0, 0, 1) = e$. This condition compactifies 
space to 5 3 , hence the field can be extended to the map : S 3 —> S 2 . Such maps belong 
to an element of tts(S 2 ) = Z, which are indexed by an integer Q. This integer is referred 
to as the topological charge. The topological charge can be found in one of two ways. 
Firstly, the tensor F^ defines a 2-form on the 3-sphere, which because H 2 (S 3 ) = is a 
closed and exact form. Therefore we can define a 1-form, A ) such that F^ = djA^ — dj-Aj 
and then Q can be found as the integral 

Q = J ZijkAFjk- (2.2) 

Generically preimages of points on S 2 will be unions of disjoint loops in M 3 . It has been 
shown [9] that Q is equal to the linking number of the loops created as the preimages of 
two distinct points. 



It can be easily seen that (2.1) is positive definite for k > and f3 < 0, this is the 



coupling region mentioned earlier and studied in [6j. It has also been shown JS] that (2.1) 



is positive definite for k < 0,/3 < with the constraint that 2/3 /k > 1. To see this we 



first re- write (2.1) as 



E = jf a 'di<f>-^ (-[FyF** - l -{d i( f> ■ d^f] + - l){di4> ■ ^0) 2 ) | d 3 x, 

(2.3) 

and, as in [TO] , we define the 3x3 matrix D ab = g l i \di<fi a )(dj<fi b ). This matrix has two 
non-zero eigenvalues which we choose to be Ai and A 2 . Then we can instantly see that 
di4> - di4> = Ai + A 2 and FijF 1 ! = 2AiA 2 , and we discover the identity 

FijF ij - \{di4> ■ d i( f>) 2 = -hx, - A 2 ) 2 . 



Therefore (2.3) is positive definite if 

k<0, (3 < 0, 2/3/« > 1. 
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This is the coupling constant region of interest in this paper. The interest in this 
region is because k and /3 have the same sign, which is needed to agree with the effective 
Lagrangian in [6j. Also, it has been shown [11] that if one restricts to the sector 2/3 /k = 1 
and introduces a constraint the model becomes integrable and exact vortex solutions can 
be constructed. 

Similar to the Vakulenko and Kapitanski argument [12] for the Skyrme-Faddeev model, 
this model also has a lower energy bound [8]. To find this bound we first use the inequality 
2FijF l i < (di4> • di4>) 2 [10J (this can be easily found from the above eigenvalues) and re- 



write the static-energy functional (2.1) as the inequality 



E > ^ J d 3 *^2]p^ - _ i)|( J d'xF^)^ (2.4) 

d 3 x^/F~^) 2 ( [ d 3 x(FijF ij ) \ 2 . (2.5) 




Now using the Sobolev-type inequality 

C (^J d'x^Ff^j d 3 xF^ > 8 ( JL J d'xe^AF^ ' , (2.6) 

where C is a universal constant and using Ward's argument [10] that C = 1/64v / 2tt 4 is 
optimal gives the bound 

E > 6^yif(^-l)|Q 3/4 . (2.7) 

Choosing a point on S 2 as the boundary value defines a unique antipodal point and 
we define the preimage of this point as the location curve. In this case the location curve 
is the curve 03 = —1 and is represented by the blue curve in all the images. The red 
curve represents the preimage of another point near to 3 = — 1 and is used to show the 
linking. 

In [8] the authors used a toroidal ansatz to reduce the (3 + l)-dimensional problem to a 
euclidean 2-dimensional problem. This greatly reduced the computational power needed 
to find minimum-energy soliton solutions. They proposed a number of minimum-energy 
solutions for (3e 2 = 1.1, but because of the ansatz these solutions where constrained to 
be planar. 

Our research is to build on this by finding minimum-energy solutions in full 3-dimensional 
space, with no constraint or ansatz. The motivation for this is because it is well known 
that most of the topological charge specific minimum-energy solutions of the Skyrme- 
Faddeev model are not planar. In fact only the topological charges one, two and four are 
planar. Hence we expected to identify lower energy solutions for the topological charge 
three and find new solutions for topological charges greater than four. Also, in the usual- 
Skyrme-Faddeev model the minimum-energy topological charge four can be understood 
as two topological charge two Hopfions on top of each other [3] , but in this model it has 
been found that the minimum-energy solution is a single topological charge four ring. 
Hence we expect to find new soliton solutions for larger topological charges. 
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3 Knotted ansatz 



To confidently identify minimum- energy solutions, we are forced to try a number of initial 
configurations with the correct topology. We can generate configurations with the correct 
topology using a technique described in [2]. We shall briefly outline this technique here. 
Firstly, we identify space with S 3 via the degree-one spherically-equi variant map, 

(Z 1 ,Z ) = ( (x 1 +ix 2 )^^,cosf + i S ^^x 3 ) , (3.8) 



where f(r) is a monotonically decreasing profile function, with the boundary conditions 

/(0)=7T,/(00)=0. 

This identifies the space as an S 3 C C 2 . We can then identify the stereographic projection, 
W, of (f) with the rational-map, p(Z , Zi)/q(Z , Z\), of two polynomials p, q. 

If we choose the polynomials such that the number of preimages of a point is Q Q we 
then have a topological charge Q configuration [2]. We can choose how the soliton 'knots' 
by choosing q(z) = to be a curve which describes an appropriate knot. The simplest 
example is the axially-symmetric planar ring which is found from the rational-map 

™=%- ( 3 - 10 ) 

This produces a field configuration where 3 = — 1 on a planar ring and has topological 
charge Q = nm. Throughout the text we categorise loop-like field configurations as A n ^ m . 
We can also create knotted field configurations with the rational map 

yet yfi 

W= ^L\ (3.11) 
Z$ + Z b y J 

This rational-map has charge Q = ab + /3a and if a and b are coprime positive integers 
the location (0 3 = —1) describes an (a, 6)-torus knot. We denote such a configuration 
as JC a ^. As well as knotted and ring configurations there are also linked solutions. We 
denote a configuration of a charge d ring linking with a charge / ring 7 times as C^'J. 
It is quite simple to see that such a configuration has charge Q = 27 + d + /. We can 
create such field configurations with a reducible rational-map, for example we can create 
a Q = 5 linked Hopfion with the map, 

W =2iZ^) + W (312) 

4 Numerical solutions 

In this section we present the numerical technique used to find minimum-energy solu- 
tions, and we discuss the minimum-energy solutions for two choices of k and /3. To 



x This is a definition of topological degree. 
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find minimum-energy solutions of the extended-Skyrme-Faddeev model we performed a 
numerical full-field minimisation. This was done using a lattice of 180 3 points with lat- 
tice spacing Ax = 0.08, this size was found to be large enough to contain the soliton 
with enough accuracy. All of the initial conditions were created using the appropriate 
rational-map, where f(r) was manually chosen to minimise the initial energy. This is 
a significantly faster technique than performing computationally intensive Hopfion colli- 
sions. 

4.1 Close to the scale-invariant limit 

In this section we chose k—— 2,/3 = — 1.1, this allows us to directly compare our results 
with those in [8j . Our numerical solutions are presented in figure [j] and the energies are 
shown in table [Q 




Figure 1: Soliton solutions of the extended-Skyrme-Faddeev model. The blue surfaces 
represent the location curve and the red is the second linking curve. Top left is charge 
one, top right is charge two, bottom left is charge three and bottom right is charge four. 

Figure [j] shows the buckled Q = 3 that we expected to discover, which could not 
be discovered in [8j and has lower energy. Also we have found a similar A^i minimum- 
energy solution to the one found in [8j, which is different to the solution found in 
the Skyrme-Faddeev model [3]. The minimum-energy solutions for Q = 1, 2, 4 slightly 
differ the the axial- ansatz solutions in [8j , the difference might be attributed to numerical 
accuracy and the extra degrees of freedom of our three-dimensional minimisations. 
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Charge 


E 


Q = l 


232.482 


Q = 2 


387.432 


Q = 3 


573.349 


Q = 4 


758.276 



Table 1: E is the minimum-energy solutions found using full three-dimensional minimi- 
sation. 

4.2 Away from the scale- invariant limit 

A feature of knotted Hopfions is that as the topological charge increases the space of 
saddle point energy solutions increases. This makes finding the minimum-energy solutions 
more involved. An important issue is that as the charge increases there are non-planar 
knotted solutions which self intersect. It was found that as the soliton tubes dynamically 
intersected they became unstable for k = — 2,/3 = —1.1. The actual cause for this 
instability is not yet understood, but it is believed to be related to the fact that at 
hi = —2,(3 = —1 the model becomes integral and has scale-invariance. Potentially the 
numerical grid is not exact enough to be stable close to this limit. Hence to proceed we 
arbitrarily fixed n = — 2,/3 = —3/2, this was found to be stable enough for the knots to 
self intersect during minimisation. This is the value of the coupling constants used in the 
remainder of this text. 



Charge 


E 


E / E bonnd 


Initial — > final configuration 


Q = l 


552.178 


1.236 


^1,1 A hl 


Q = 2 


903.544 


1.203 


-4 2 ,i ->• -4 2 ,i 


Q = 3 


1267.24 


1.245 


-4-3,1 ->• -4 3 ,i 


Q = 4 


1606.32 
1687.9 


1.272 


-4.4,1 ->• -4 4 ,i 

L{i, ^2,2 - >■ -42,2 


Q = 5 


1958.73 


1.130 


-45,1, ^3,2, L{ 2 -45,1 


Q = Q 


2263.04 


1.322 


-46,1 j ^3,2 — > L 2 'o 


Q = 7 


2506.66 
2528.81 
2647.1 


1.304 


r 2,3 r 2,3 
^-1,1 — > -^1,1 

^3,2 - >■ ^3,2 

-4 7 ,i ->• -4 7 ,i 



Table 2: The minimum-energy found using full three-dimensional minimisation for k — 
— 2,/3 = —3/2 and the initial and final forms. 

The minimum-energy solutions we have found for Q = 1, 2, 3 are similar to the Q = 
1,2,3 minimum-energy solutions of the Skyrme-Faddeev model. As shown in figure [2d] 
the Q = 4 minimum-energy solution is a larger ring, this is different to the A 2 ^ Q = 4 
minimum-energy Skyrme-Faddeev solution, also for f3 = — 3/2, k = —2 the ring is buckled. 
The f3 — — 1.1, ft = — 2Q = 4 is not very buckled. We have also found linked (5 = 4 
solutions which are energetic local minimum or saddle point solutions. From this we 
can conject that the extended-Skyrme-Faddeev model energetically prefers large vortex 
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(a) Q = 1 



(b) g = 2 a 2( i 



(c) Q = 3 A 3 ,i 






(d) Q = 4 A 4 ,i -> A 4 ,i 



(e) Q 



4 Lft -> A 2 ,2 



(f) Q = 4 A 2)2 -»■ -42,2 




(g) Q = 5 K 3 ,2 -> 




(j)Q 



7 Ln o — ^ -Zvo' 



J 2,3 



-2,3 





(h) Q = 5 L\l A 5 ,l 



(i) Q = 6 A 6 ,x L\\ 





(k) Q = 7 /C 3 , 2 /C 3>2 



(1) Q = 7 A 7jl 



Figure 2: Hopfions of the extended-Skyrme-Faddeev model for k = — 2,/5 
each Q the upper most left is the minimum-energy solution we found 



-47,1 

-3/2. For 



tube configurations and not composite configurations. Again the Q = 5 minimum-energy 
solution is a large ring (figure 2gf2h ), where we have found that a knotted (^3,2) and a 
linked (L\\) initial configurations have deformed to a minimum-energy ring-like (A^i) 
solution. Also, as shown in table [2j the Q = 5 minimum-energy is much closer to the 
lower bound than either the Q = 4 or Q = 6. The Q = 7 minimum-energy we find is a 
linked Q = 3 and Q = 2 (figure |2jj) , this solution has E = 2506.66 and the solution 



has E = 2528.81 (figure 2k). These energies are too close to be able to confidently identify 
the minimum-energy solution. 
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5 Concluding remarks 

We have presented new minimum-energy solution candidates for the extended-Skyrme- 
Faddeev model. We have found that for 1 < Q < 7 the minimum-energy solutions 
are large single-core Hopfion loops. This is notably different from the minimum-energy 
solutions of the Skyrme-Faddeev model. This difference is not surprising as our coupling 
is very different to the Skyrme-Faddeev model, which is known to possess linked and 
knotted solutions for Q < 7. From the minimum-energy solutions found it seems that 
the extended-Skyrme-Faddeev model energetically prefers large buckled ring-like field 
configurations. It has recently been shown that the minimum-energy solutions of the 
Skyrme-Faddeev model can be approximated by elastic theory [15] . Here the authors 
understood the knotting and linking form of minimum-energy solutions of the Skyrme- 
Faddeev model with a minimal elastic-rod model. They showed that for Q > 2 the 
minimum-energy solutions are buckled rings. Here we expect the same elastic-rod model 
to replicate the minimum-energy solutions, but with a different coupling parameter and 
to predict the buckling. This would be very interesting to study. 
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